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ABSTRACT 

We present an analysis of the luminosities and equivalent widths of the 284 z < 0.56 [O II]-emitting 
galaxies found in the 169 arcmin^ pilot survey for the Hobby-Eberly Telescope Dark Energy Exper- 
iment (HETDEX). By combining emission-line fluxes obtained from the Mitchell spectrograph on 
the McDonald 2.7-m telescope with deep broadband photometry from archival data, we derive each 
galaxy's de-reddened [O II] A3727 luminosity and calculate its total star formation rate. We show 
that over the last ^ 5 Gyr of cosmic time there has been substantial evolution in the [O II] emission- 
line luminosity function, with L* decreasing by ~ 0.6 ± 0.2 dex in the observed function, and by 
~ 0.9 ± 0.2 dex in the de-reddened relation. Accompanying this decline is a significant shift in the 
distribution of [O II] equivalent widths, with the fraction of high equivalent-width emitters declining 
dramatically with time. Overall, the data imply that the relative intensity of star formation within 
galaxies has decreased over the past '^ 5 Gyr, and that the star formation rate density of the universe 
has declined by a factor of ~ 2.5 between z ~ 0.5 and z ~ 0. These observations represent the 
first [O II]-based star formation rate density measurements in this redshift range, and foreshadow the 
advancements which will be generated by the main HETDEX survey. 

Subject headings: galaxies: formation — galaxies: evolution — galaxies: luminosity function - cos- 
mology: observations 



1. INTRODUCTION 

The evolution in the cosmic star formation rate den- 
sity (SFRD) is an important observational constraint for 
all the cu rrent models of galaxy formation and evolu- 
tion (e.g.. lSomerville et al.ll2012i) . It is generally agreed 
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that the SFRD reached its peak around z ^ 2, and since 
then has declined b y roughly an orde r of magnitude (e.g., 
Madau et allll996l: iLillv et al.lll996l: IHopkins fc BeacomI 
i2006i : iBouwens et al.ll2010l ). However, the details of this 
evolution remain poorly understood, in part because of 
the systematic uncertainties associated with the patch- 
work of various star formation rate (SFR) indicators. 

For example, in the nearby universe, the strength of 
the Ha emission line is one of the most direct and robust 
measures of ongoing star formation, since it is powered by 
the photoionization produced by massive (M > IOMq), 
young (t < 20 Myr) stars (Kennicutt 1991). Unfortu- 
nately, at distances greater than z ~ 0.4, this line red- 



shifts out of the optical and into the near-IR portion 
of the spec trum, where it is much more difficult to ob- 
serve (e.g.. iGlazebrook et al.l [200^ . Consequently, by 
z r^ I, the [O II] A3727 emission line, which is produced 
by coUisional ex citation, has taken the place of Ha (e.g., 
iHippelein eTaLi 2003: Lv et al. 2007 ). and by ^ > 2, the 
preferred indicator of star formation is the rest- frame 
ultra violet continuum (e.g., iPettini et al.ll200ll: iHopkinsI 
l2004f ). Like the emission-line indicators, the rest- frame 
UV directly traces the flux of young, massive stars, but 
in this case, the timescale over which the star forma- 
tion is measured is ^ 100 Myr, and the result is much 
more sensitive to internal extinction. Mid-IR measure- 
ments, which reflect the re-radiation of stellar emission 
by warm dust, integrate s tar formation over an even 
longer (^^ 1 Gyr) timescale (jKelson fc Holdenll2010[ ). fur- 
ther complicating the interpretation of the low-z epoch, 
when the SFRD of the universe is evolving rapidly. 

In practice, a large number of additional star formation 
rate indicators are used both at high and low redshift, 
and these span the entire electromagnetic spectr um, from 
1.4 G Hz radio and far-IR bands to the X-ray (seelHo_p.kiri^ 
12004 and references therein) . One of the most versatile 
of these is the flux from the [O II] doublet at 3727 A. 
As a SFR indicator, [O II] has two significant advan- 
tages over measurements in the near and far UV: it is 
easier to detect at low and intermediate (z < 2) redshifts 
and it is more sensitive to the instantaneous star forma- 
tion rate, rather than the ~ 100 Myr time-averaged rate. 
Unfortunately, the path from [O II] emission to star for- 
mation rate is much less direct than for either Ha or the 
near-UV, as the emission depends on parameters such as 
the galactic oxygen abundance and the ionization state 
of the gas (Kcnnic utti Il998l ). Thus, the method must 
be calibrated er npirically via comparisons to other SFR 
indicators fe.g.. JJansen et al.ll200lt iKewlev et al.l 1200^ 
iMoustakas et al.ll2006D . 

Over the past 15 years, there have been several [O II]- 
based measurements of SFRDs, mostly at z ^ and 
z ~ 1. Many of these studies derive from spectroscopic 
observations of magnitude-limited samples of galaxies, 
and are thus severely incomplete at the faint end of 
the emission-line luminosity function ()Hogg et al.l 119981 : 
iHammer et"al]|1997l: [Zhu et al.ll2009D . Others are based 
on narrow-band surveys of selected high redshift {z > 
0.8) epochs, and are thus only sensitiv e to high equiva- 
lent w idth objects fe.g.. ILv et al...2007t iTakahashi et al.l 
|2007() . Still others employ near-IR Fabry-Perot observa- 
tions and are restricted to the very brightest (^ 1 dex) 
[O II] emitters at z ~ 1 ()Hippelein et al.ll2003[) . Con- 
sequently, while the local (z < 0. 2) [O II] luminosity 
function is reasonab ly well defined (jGallego et al. 199g 
IGilbank et al.|[2"010D . measurements of the evolution of 
[O II] A3 72 7 emission over the last ~ 5 Gyr of cosmic 
time are rather limited. In fact, the only studies that 
attempt to sample the full-range [O II] over these inter- 
mediate redshifts are the Hubble Space Telescope grism 
surveys of lTeplitz etHI (|2003|) and JPirzkal et all (|2012[ ). 
Unfortunately, these z > 0.5 observations are limited by 
a rather high equivalent width limit (> 35 A), and a lack 
of information concerning the galaxies' internal extinc- 
tion. 

There is, however, one other dataset that can be ex- 



ploited to measure the evolution of [O II] emission in the 
nearby universe. The pilot survey for the Hobby-Eberly 
Telescope Dark Energy Experiment (HETDEX) ob- 
tained blind integral field spectroscopy over 169 arcmin^ 
of sky and measured [O II] A37 27 fluxes for 2 84 galax- 
ies between < z < 0.56 (Adams et al.ll201lD . Only a 
few [O II] surveys cover this redshift range, and most 
are defined by magnitude-limited samples of objects. In 
contrast, the HETDEX pilot survey's dataset is emission- 
line flux limited, rather than continuum flux limited, and 
is thus complete to much lower star formation rates than 
surveys which choose their targets via broadband mea- 
surements. The observations therefore sample a section 
of parameter space rarely probed by other investigations, 
and is a unique resource for studying the relationship 
between [O II] and other star formation rate indicators. 
Moreover, because the HETDEX pilot survey was con- 
ducted in fields with large amounts of ancillary data, we 
also have access to each galaxy's spectral energy distribu- 
tion (SED), enabling a full examination of its photomet- 
ric properties and internal reddening. Finally, because 
the survey's redshift range covers the last ~ 5 Gyr of 
cosmic time, the database can be used to trace the late- 
time evolution of [O II], and provide a comparison for 
SFRD measurements obtained at higher redshift. 

In this paper, we present an analysis of the [O II] A3727 
emission-line galaxies found in the HETDEX pilot sur- 
vey. In §2, we summarize the observations and the se- 
lection method which produced the catalog of 284 [O II] 
line emitters. In §3, we present the evolution of the [O II] 
equivalent width distribution with redshift, and use this 
reddening-independent quantity to show how the rela- 
tive importance of star formation has declined over the 
last ^ 5 Gyr. In §4, we examine the internal extinction 
of the [O II] emitters, and compare the star formation 
rates derived from the measurements of the [O II] flux to 
those determined via the rest-frame UV continuum. In 
§5, we present both observed- and extinction-corrected 
[O II] luminosity functions and explore the evolution of 
the cosmic SFRD out to z ~ 0.5. We conclude by dis- 
cussing the implications of our work for future surveys of 
star-forming galaxies. 

Throughout this work, we adopt a ACDM cosmology 
with ^m = 0.3, f^A = 0.7 and Hq = 70 km s"^ Mpc'^ 

2. THE HETDEX PILOT SURVEY DATA 

The data for our analysis were taken as part 
of the HETDEX pilot survey, a blind integral-field 
spectroscopic study of four areas of sky containing 
a wea lth of ancillary m ulti-waveleng th data: COS- 
MOS (IScoville et ani2007h. GOO DS-N (IGiavalisco et al.l 
[200l. MUNICS-S2 (|Drorv et al.l[200l . and XMM-LSS 
( Pierre et al.l I2004D . A complete description of this 
pathfinding survey and its data products is given in 
I Adams et al.l (|2011D . Briefly, a square 246-fiber array 
was mounted on the Harlan J. Smith 2.7-m telescope at 
McDonald Observatory, and coupled to the George and 
Cynthia Mitchell Spectrograph, a proto-type of the Visi- 
ble Integral-field Replicablc Unit S pectrograph (VIRUS) 
designed for HETDEX (Hil l etTaLl l2008bl This instru- 
ment was set to disperse the wavelength range between 
3600 A and 5800 A at -- 5 A resolution and 1.1 A pixel" ^ . 
Each fiber subtended 4'.' 2 on the sky; consequently, the 
entire fiber array covered an area of 1.'7 x 1.'7 with a 1/3 
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fill factor. To fill in the fiber gaps and improve the sur- 
vey's spatial resolution, each field was observed using a 
6-pointing dither pattern, generating 6 x 246 = 1,476 
separate spectra. In total, the pilot survey covered 
169 arcmin^ of blank sky (71.6 arcmin^ in COSMOS, 
35.5 arcmin^ in GOODS-N, 49.9 arcmin^ in MUNICS- 
S2, and 12.3 arcmin^ in XMM-LSS) with 1 hour of ex- 
posure time at each dither position. Fifty percent of 
the pointings reached a 5 cr monochromatic flux limit of 
6.7 X 10-1^ ergs cm^^ g-i ^t 5000 A, and 90% reached 
1.0 X 10~^^ ergs cm~^ s^^. Simulations demonstrate that 
for objects with equivalent widths greater than ^^ 5 A, 
the recovery fraction of emission lines integrated across 
all continuum brightnesses is greater than 95%, and even 
for equivalent widths as small as 1 A, the completeness 
fraction is better than 90%. 

The goal of the HETDEX pilot survey was to test 
the equipment, analysis procedures, and ancillar y data 
requir ements for the main HETDEX survey (Hi ll et al.l 
l2008al) by performing an unbiased search for emission- 
line objects over a wide range of redshifts. To do 
this, each sky-subtracted fiber spectrum was first flux- 
calibrated and then searched for the presence of emis- 
sion lines. When a candidate emission feature was found, 
the COSMOS, GOODS-N, MUNICS-S2, and XMM-LSS 
imaging data were used to identify the line's photometric 
counterpart. The combination of the spectroscopic line 
flux and the photometrically determined continuum flux 
density were then used to calculate the line's equivalent 
width. The identity of the emission line and the object's 
most-likely redshift were then determined using either 1) 
the presence of multiple emission lines in the spectrum, 
2) the source's photometric spectral-energy distribution 
(as determined from the multi- wavelength images), or 3) 
the emission-line's rest-frame equivalent width. The lat- 
ter criterion is particularly useful for discriminating Lya 
from local [O II] A3727, as single-line sources with rest- 
frame equivalent wid ths greater than 20 A are almost 
always Lya emitters (jGronwall et al.|[2007D . 

In total, the HETDEX pilot survey identified 397 
emission-line sou rces, including 104 Lya emitters 
(jBlanc et al.l 1201 1[ ). and 284 [O II] emission line galax- 
ies. The latter have redshifts between 0.078 < z < 0.563, 
[O II] emission- line luminosities between 39.8 < logi < 
42.2 (ergs cm~^ s~^)j and [O II] rest-frame equivalent 
widths as great as 77 A. 

3. STAR FORMATION VERSUS EMISSION FROM AGN 

Before examining the issue of star formation in our 
emission-line galaxies, we first must consider the ef- 
fect that active galactic nuclei have on our [O II] 
fluxes. The pilot survey spectra do not have the wave- 
length coverage to employ traditional line-ratio diagnos- 
tics s uch as [N II] /Ha or [ O III]/H/? for AGN detection 
(e.g.. iBaldwin et all [l98lt iVeilleux fc Osterbrockl [19871: 
iKewlev et al.l 120011 ^ However, because the fields cho- 
sen for study have been studied with the Chandra and 
XMM observatories, we can identify probable AGN con- 
tami nants via X-r a y lurn inosity. 

In lAdams et al.l (|2011[ ) , X-ray counterparts were iden- 
tified for 30 of the HETDEX pilot survey's [O II] emit- 
ters. However, many of these objects are in the GOODS- 
N region, where the Chandra images are deep enough to 



detec t the X-rays associated with normal sta r formation 
(e.g., iPersicet al.|[20p: iLehmer et aL|[20lol ). If we ex- 
clude the [O II] emitters with luminosities fainter than 
10'*° ergs s~^ in the 2 to 8 keV X-ray band, we are left 
with 10 probable AGN: 5 in GOODS-N (which has the 
deepest X-ray coverage), and 5 in the other three fields. 
This number, which represents ^ 3% of the sample, is 
slightly greater tha n the AGN contamination rate of 1 
to 2% estimated bv iZhu et al.l ([2009) in their analysis of 
the line-ratios and X-ray emission of z '--^ 1 [O II] emitters 
in the DEEP2 survey. Since the X-ray luminosity limits 
used here are more tha n an order of ma gnitude fainter 
than those considered bv lZhu et al.l (|2009( ). we judge this 
agreement to be acceptable. 

Of course, just because an object emits X-rays, it does 
not necessarily follow that the system's integrated [O II] 
emi s sion i s dominated by flux from the central engine. 
iHol (|2005[ ) has shown that the physical conditions in 
the narrow-line regions of AGN disfavor the creation of 
strong [O II]. Indeed, an examination of the X-ray bright 
objects in our survey demonstrates that these likely AGN 
are generally not amongst the brightest [O II] sources. 
Consequently, their existence does not substantially al- 
ter the observed [O II] luminosity function or equivalent 
width distribution, nor do they effect the main conclu- 
sions of this paper. We therefore exclude these AGN 
candidates from our analysis. 

4. THE EQUIVALENT WIDTH DISTRIBUTION 

One quick but very coarse way of examining the rela- 
tive importance of star formation as a function of time is 
through the use of the distribution of [O II] A3 72 7 rest- 
frame equivalent widths. The [O II] flux from a galaxy 
generally measures the amount of star formation in the 
very recent past, i.e., t < 20 Myr (|Kennicuttlll998[ ). In 
contrast, the continuum underlying 3727 A in large part 
reflects st ar formation over a much longer (~ 1 Gyr) 
timescale (jGilbank et al.l 120101 ). Thus, by forming the 
ratio of these two quantities, it is possible to examine the 
relative importance of the current burst of star formation 
in a manner that reduces the complications introduced 
by internal extinction. 

In the local neighborhood, the distribution of [O II] 
A3727 rest-frame equivalent widths peaks near '-~^ 5 A, 
and then rapidly decays, as would an exponential with 
a scale length of wq < 10 k (|Blanton fc L in 2000). 
Notably, this function has virtually no dependence on 
the brightness of the galaxy: the distribution of rest- 
frame [O II] equivalent widths for the ~ 8500 galaxies 
discovered in the Las Campanas Redshift Su rvey is al- 
most independent of -R-band luminosity (Blan ton fc LinI 
|2000( ). Similarly, there is no evidence for a correlation 
between luminosity and equivalent width in an Ha grism- 
survey of the local un i verse : if we divide the 191 galaxies 
of the iGallego et al.l (119961) survey in half and per form 
an Anderson-Darling test (jScholz fc StephensI [l987l ) . we 
find no statistical difference between the distributions of 
[O II] equivalent widths for the high- and low-line lumi- 
nosity samples. 

Figure [T] displays the rest-frame [O II] equivalent width 
distribution for the HETDEX pilot survey's sample of 
emission-line galaxies, as a function of redshift and [O II] 
luminosity. Specifically, the figure sub-divides our sam- 
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Figure 1. The distribution of [O II] A3727 rest-frame equivalent 
widths for the pilot survey galaxies, plotted as a function of [O II] 
luminosity (j/-axis) and redshift (x-axis) . The location of each his- 
togram illustrates the redshifts and (approximate) luminosities of 
the galaxies contained within it. The dotted red line shows the 
80% detection limit, i.e., 80% of the survey frames have 5cr detec- 
tion limits brighter than this threshold. The [O I] airglow feature 
at 5577 A is easily seen. A total of 259 galaxies lie above the 80% 
threshold. The curves display how the best-fit exponential changes 
with redshift. Since the distribution of equivalent widths turns 
over at small values, objects with rest-frame equivalent widths less 
than 5 A have been excluded from the calculation. The absence of 
high-luminosity [O II] emitters at low redshift is partially a volume 
effect. Note the distribution of equivalent widths is insensitive to 
[O II] luminosity, but is strongly dependent on redshift. 

pie of [O II] emitting galaxies into four redshift bins, and 
histograms the rest-frame equivalent widths for galax- 
ies in the top-half and bottom-half of the emission-line 
luminosity function. The location of each histogram il- 
lustrates the redshifts and (approximate) luminosities of 
the galaxies contained within it. The dotted line in the 
figure illustrates the data's "80% completeness limit," 
i.e., 80% of the fields observed in the survey have 5cr de- 
tection limits brighter than this value. Only one out of 
the 42 galaxies dX z < 0.2 has a line-luminosity brighter 
than logL = 41.1; this is partially a volume effect, as 
the probability of finding a high-luminosity target in the 
nearby universe is small. Conversely, at z > 0.45, low- 
luminosity [O II] emitters fall below the flux limit of the 
survey. Nevertheless, despite this selection effect, it is 
clear that any dependence of equivalent width on [O II] 
luminosity is minimal. While we cannot exclude the pos- 
sibility of a small luminosity dependence in the equiva- 
lent width distribution, an Anderson-Darling test finds 
no significant difference between the high-luminosity and 
low-luminosity samples, even at the Icr level. 

Conversely, the redshift dependence of the rest-frame 
[O II] equivalent width distribution is quite strong, with 
the number of high-equivalent width objects increas- 
ing rapidly with redshift. At z < 0.2, the median 
rest-frame [O II] equivalent width is ~ 10 A, a value 
only slightly larger than the 7.5 A median found lo- 
cally from the m agnitude-limited Las Campanas Survey 
(jBlanton fc Linli20003 . By z > 0.45, however, the median 
equivalent width has increased to ^ 21 A. Qualitatively, 
this shift is similar to that seen in the magnitude-sel ected 
samples of galaxies o bserved by iHogg et all ()1998D and 
I Hammer et al.l (|1997f ). although in these studies, most of 
the increase occurs at z > 0.5. The slitless HST mea- 
surements of iTeplitz et al.l ( |2003,) imply a much larger 
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Figure 2. Relative likelihoods that the rest-frame [O II] equiv- 
alent widths plotted in Figure [T] are drawn from an exponential 
distribution with scale length Wo- Objects with equivalent widths 
less than 5 A have been excluded from the fit. In the local uni- 
verse, rest-frame [O II] equivalent widths e-fold with a scale length 
of ~ 8 A; by z ~ 0.5, this scale length has increased to ~ 22 A. 

equivalent width scale length (^ 50 A) with no signifi- 
cant evolution between 0.5 < z < 1.2. However, since 
these parallel STIS observations were only sensitive to 
objects with EWobs > 35 A, our new measurements are 
n ot necessarily in confli ct with their result. 

iBlanton fc LinI (|2000r ) have shown that for a local, 
magnitude-limited sample of galaxies, the distribution 
of [O II] equivalent-widths is well-fit with a log-normal 
function with a mean of (EWq) '^ 10 A and dispersion 
of cr ^ 0.77 A. When we fit our emission-line selected 
data in this way, we obtain a slightly larger value for 
the mean of the distribution, (EWq) ~ 15 A for all but 
the highest redshift bin, and a standard deviation that 
gradually increases from tT~0.5Aat2:<0.25tO'^lA 
at z '^ 0.5. Such a parameterization does reproduce the 
turnover seen at small equivalent widths, while maintain- 
ing the high-equivalent width tail of the distribution. It 
is, however, not the most intuitive of functions, especially 
for visualizing the evolution of the strongest [O II] line 
emitters. Thus, to quantify the evolution of the high- 
equivalent width objects displayed in Figure [TJ we chose 
to use a simple exponential law. We excluded those ob- 
jects with equivalent widths less than 5 A, combined 
the low and high-luminosity samples at each redshift, 
and computed the relative likelihood that the observed 
data were drawn from a model exponential with scale 
length. Wo- In order to account for the heteroskedastic 
nature of the dataset, we performed this experiment 1000 
times, using different realizations based on the quoted 
uncertainties of each measurement. As the probability 
curves of Figure [5] demonstrate, the scale length of the 
rest-frame [O II] equivalent width distribution increases 
smoothly with redshift, going from Wq = 8.0 ± 1.6 A 
at z < 0.2 to Wo = 21.5 ± 3.3 A at z > 0.45. The 
simple interpretation of this trend is that over the past 
'^ 5 Gyr, the relative intensities of individual starbursts 
has decreased linearly with time. 

5. INTERNAL EXTINCTION WITHIN THE SAMPLE 

Before we can translate the HETDEX [O II] lumi- 
nosities into star formation rates, we must first examine 
the effect that internal extinction has on each galaxy's 
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Figure 3. SED-based stellar reddening values for the 259 VIRUS- 
P [O II] emitting galaxies above tlie 80% completeness limit, as a 
function of function of galaxy absolute magnitude. The grey area 
approximates the locus of points found by Moustakas et al.i (|200Q) 
using the Balmer decrements of a local sample of galaxies. The 
agreement between the two samples is excellent, and there is no 
evidence for a redshift dependence in the extinction. The crosses 
represent seven probable AGN. 



emergent emission-line flux. Often times, this is done 
by applying a mean reddening corr ection to the pop- 
ulation of galaxies as a whole (e.g., iFuiita et al.l 120031 : 
ILv et al - 2007): in other cases, a mean extinction is de- 
fined as a function of galaxy abso lute m agnitude (e.g., 
iMoustakas et al.l 120061 : iZhu et all l2009( ). However, a 
more robust way of handling the problem is to de- 
redden each galaxy individually using information from 
the galaxy itself. Although the HETDEX spectra lack 
sufficient wavelength coverage to constrain internal ex- 
tinction via the Balmer decrement, all of the survey fields 
have been imaged in a variety of bandpasses, from the 
UV (with GALEX) through the near-IR. Thus, it is pos- 
sible to obtain an estimate of each object's internal ex- 
tinction by fitting its spectral energy distribution (SED) 
to models of galactic evolution. 

To do this, we u s ed th e population synthesis code of 
iBruzual fc CharloO ()2003[ ) to create a three-dimensional 



grid of models us ing the assumptions of solar-metallicity 
stellar evolution (jBressan et al.lll993D. alSalpeteii ()1955[ ) 
initial mass function, and a ICalzettil ()2001[ ) reddening 
law. In total, 504 models were produced, using four 
different values for the e-folding timescale of star for- 
mation (r =10, 50, 100, and 500 Myr), 21 logarfth- 
mically spaced galactic ages (1 Myr < i < 10 Gyr), 
and six uniformly spaced values for the stellar reddening 
(0 < E{B -V) < 0.5). We then found the best-fitting 
values of r, t, and E{B — V) by fitting each galaxy's ob- 
served multicolor photometry to the computed colors of 
the grid via a x^ minimization code. Obviously, single- 
component population synthesis models such as these 
are somewhat limited, as they cover only a small part 
of galaxy parameter space. Nevertheless, the models do 
provide some guidance as to a galaxy's internal extinc- 
tion, and therefore represent an important improvement 
in the study of the evolution of star formation. 

Figure [3] shows the SED-based stellar reddenings for 
our sample of [O II] emitters, as a function of galactic 
absolute magnitude and redshift. For the figure, we have 
adopted the extinction curve of .Calzc tti (.2001i ). along 
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which ICalzettil (|2001l ) inferred from UV, far-IR, and 
Balmer-line observations of ^ 50 starburst galaxies in the 
local neighborhood. The figure confirms the well-known 
correlation betwe en a galaxy's absolute m a gnitude and 
extinction (e.g., 'Wang fc Heckmaiil 119961 : iTullv et al.l 
[1998: Jansen et al. 2001). Indee d, the use of equat i on ([T| ) 
recovers the relation observed bv lMoustakas et aD ()2006f ) 
for galaxies in the local neighborhood. Also, this is no 
clear evidence for a shift in the relation with redshift, 
although the size of the dispersion {(Te^b-v) ~ 0-2), the 
limited number of objects, and the dependence of sur- 
vey volume and completeness on redshift all conspire to 
make measurements of the evolution of extinction diffi- 
cuft. The full sample of ~ 10^ HE TDEX galaxies w ill 
greatly improve upon this situation (|Hill et al.ll2008aD . 

Another test of the robustness of our reddening esti- 
mates is to compare the de-reddened [O II] A3727 lu- 
minosities with the galaxies' de-reddened near-UV mag- 
nitudes, as recorde d in the GR6 catalog of GALEX 
([Martin et al.l 120051 ). Both datasets probe the presence 
of hot, young stars: [O II] through the indirect detection 
of ionizing radiation, and the UV via the direct measure- 
ment of the st ellar continuum of hot stars. However, ac- 
cording to the ICalzet ti (2001) extinction law, the extinc- 
tion that effects [O II] A3727 emission line is ~ 3.5 times 
greater than that which extinguishes a galaxy's UV con- 
tinuum. Consequently, any systematic error in our red- 
dening estimates should be easily visible in a comparison 
of the inferred star formation rates. 

To perform this test, we u sed the UV conti nuum star 
formation rate calibration of lKennicuttI ([19981 ) 



SFR(UV)(Mq yr^i) = 1.4 x 10"^*^ L^ (ergs s^^Hz 



-lxJ„-l^ 



(2) 



and the empirica l [O II] A3727 calibration of 
IKewlev et aT] ([200l 

SFR([0 II])(M0yr-i) = Q.bSxlQ-'^'^ L{[0 II]) (ergs s"^) 

(3) 
Both of these calibrations assume a iSalpeteii (|1955l ) ini- 
tial mass function truncated at IOOMq and present-day 
solar metallicity. (Because of our limited wavelength cov- 
erage, no i?23 abundance correction could be applied to 
the [O II] lumin osities.) If our r eddening estimates are ro- 
bust, and if the ICalzettil ([200 If ) extinction law holds, the 
ratio of these two star formation rate indicators should 
be one. 

Figure |4| performs this comparison. In the figure, an 
SED-based /c-correction has been applied to the GALEX 
near-UV magnitudes, and both quantitie s have be e n cor - 
rected for internal extinction using the ICalzettil ([2001f l 
reddening law. From the figure, it is clear that the 
two star formation indicators are strongly correlated. 
The best-fit slope of the relation, as determined by the 
Buc kley- James and E M algorithms for censored data 
(see llsobe et all Il986f ). is just slightly less than one 
(0.95ib0.04), \ yith t he offset in the direction predicted by 
IKewlev et al.l (|2004f ) for high star formation rate systems. 
Interestingly, the scatter about this line, (cr ~ 0.3 dex) is 
three times larger than the typical ~ 0.1 dex errors as- 
sociated with our [O II] and UV measurements, and sig- 
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Figure 4. A comparison of star formation rates derived from mea- 
sur ements of the UV continuum (using the calibration oflKennicuttJ 
119981 and [O II] flux (via the calibration of Kcwlev ct al. (2004) 
without the i?23 correction for oxygen abundance). The UV fluxes 
have been fc-corrected, and both quantities have been corrected 
for internal ex tinction via our SED-based reddening values and a 
ICalzettJ 11200 II ) reddening law. Open circles mark objects that are 
not in the GALEX GR6 catalog, and denote approximate upper 
limits. The crosses show X-ray laright objects. The scatter in the 
diagram is substantial, but any systematic error between the two 
SFR indicators is small. 

nificaiitly greater than the ~ 0.08 dex dispersion seen in 
local comparisons of Hct and [O II] star formation rates 
(e.g., iKewlev et al.ll200^ . Part of this scatter is likely 
due to our use of the mean lKewlev et aTl ()2004D calibra- 
tion for [O II] emission; without spectral measurements 
at [O III] A5007 and II/3, we cannot apply second-order 
corrections to compensate for the ionization state and 
oxygen abundance of the gas. In addition, the extinction 
that affects the forbidden [O II] emission line is not neces- 
sarily that which extinguishes the UV continuum: while 
equation ([1]) works well in the mean, there should be a 
non-negligible amount of dispersion about this relation. 
Thus, even if our SED-based extinction estimates were 
perfect, our error budget would still have this additional 
reddening term. Finally, some component of the scatter 
may be driven by the stochastic nature of star forma- 
tion. Emission-line SFR indicators, such as Ha and [O II] 
A3727, record star formation over the past ~ 20 Myr, 
while measurements of the stellar UV continuum prob e 
a timescale that is ~ 5 times longer (jKennicuttl 119981 ). 
The scatter between the two indicators must necessarily 
be larger than that seen in an IIq;-[0 II] comparison. In 
any case, the figure confirms that our estimates for inter- 
nal reddening are reasonable and should not inject any 
large systematic effects into our analysis. 

6. THE [O II] EMISSION-LINE LUMINOSITY FUNCTION 

To measure the [O II] emission-line luminosity func- 
by f ollowing the procedures described 



ti on, we bega n 

inlBlanc et al.l ( 20 111) and applied the 1 /T4nax technique 



dSchmidt 1968; IHuchra fc SargentllTOTl to the 360 sep- 



arate dithered pointings of the HET DEX pilot survey. 
As described in lAdams et al.l ()2011l) , the noise charac- 
teristics of each survey frame were recorded as a func- 
tion of wavelength. Consequently, for each [O II] emitter 
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Figure 5. The observed (left) and de-reddened (right) [O II] 
A3727 luminosity functions as a function of redshift. The black 
points in the left-han d diagram display the local luminosity func- 
tion oflGilbank et al.l (f2bl0). The dashed curves show the best-fit 
ISch echtcr (1976) function for the lowest redshift bin, while the solid 
lines show the best-fit functions at each redshift. In both plots, the 
faint-end slope of the fitted luminosity functions has been fixed at 
a = —1.2. While there is evolution in both functions, the change 
of L* with redshift is much more pronounced when reddening is 
taken into account. 



in the sample, we could calculate Knax, the co-moving 
volume of all redshifts at which the object could have 
been detected with a signal-to-noise ratio above a given 
threshold. The implied number density of galaxies in any 
absolute luminosity bin, dlogL, is then 



<^(logL) 



1 



A(logL) 



^E 



1 



K. 



(4) 



where 77 is the inverse of the completeness function, and 
the summation is performed over all galaxies with lumi- 
nosities falling within the bin. Figure [5] displays these 
derived [O II] luminosity functions (both observed and 
de-reddened) for the redshift intervals < z < 0.2, 
0.2 < z < 0.325, 0.325 < z < 0.45 and 0.45 < z < 0.56. 
In each panel, we have excluded the likely AGN (3 ob- 
jects in panels 1, 2, and 4, and 1 object in panel 3), 
binned the data into AlogL = 0.2 intervals, and drawn 
error bars based on the bin's counting statistics. 

As can be seen from the top segments of the figure, 
our z < 0.2 observed luminosity function is in very 
good agreement wi th that derived for z < 0.2 galaxies 
in SDSS Stripe 82 (IGilbank et al.l[20Tol) . In this lowest 
redshift bin, our survey volume is only ^ 2600 IVIpc'^, 
hence intrinsically rare objects which populate the ex- 
treme bright-end of the luminosity function are poorly 
represented. However, at fainter line-luminosities our 
function is reasonably well- fixed at a level t hat is ~ 40% 
higher than that found bv IGilbank et al.l ()2010[ ). This 
is consistent with the expectations, as for the galactic 
number densities under consideration, the effec ts of cos- 
mic variance should be of this order (jTrenti fc Stiavelhl 
I2OO8I) . 

At redshifts z > 0.2, there is clear evidence for evolu- 
tion, both in the observed luminosity function and, more 
dramatically in the de-reddened function. To quantif 
these changes, we fit the data of Figure [5] to a .Schechte 
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(|1976D function, i.e. 



(j}{L/L*)d{L/L*) = 0* (L/i*)" e-^/^ d[L/L*) (5) 

and examined the behavior of L* and 0* as a function 
of redshift. While o ther forms of the lumi nosity function 
are possible ~ both iZhu et all ()2009[ ) and 'Gilban k et al.l 
(jlOlO) use double power laws for their [O II] datasets - 
our sample of objects is not large enough to warrant the 
introduction of additional parameters into our analysis. 
Quantifying the evolution of the bright, intrinsically rare 
[O II] emitters requires a much larger survey volume, 
such as that which will be covered by the full HETDEX 
program. 

To fit the luminosity functions of Figure El we avoided 
the limitations associated with our finite-sized bins, and 
used a maximum-likelihood analysis. We began by com- 
puting 0'(i), the true luminosity function (in galaxies 
per cubic Megaparsec) from which the data of Figure [5] 
are drawn, modified by observational selection effects, 
such as photometric errors and incompleteness. By defi- 
nition, in a given volume of the universe, Al^, and within 
a given luminosity interval, Ai, we should expect to ob- 
serve A = (f)'{L)/S.L /S.V galaxies. From simple Poisson 
statistics, the probability of actually observing n galax- 
ies in that same interval is 



p(n\\) 



A"e~ 



(6) 



If we shrink the bin size down to zero so that the interval 
becomes a differential, then each division of the luminos- 
ity function will contain either zero or one object. The 
probability of observing any given luminosity function is 
then 



A?e-^' 



n 0! 
n ^ ^ n^^ 



n 



bins with bins witli 1 

N 

-A; 



1! 



(7) 



all bins 



1=1 



where N is the total number of galaxies observed. With 
a little manipulation, this becomes 

N 

P = eii-p{-^4>'{L)dLdv\ Wcj)' {L{)dL dV (8) 



or, in terms of relative log likelihoods 



InP 



I-Z2 /-oo N 

/ / (t>'{L)dLdV + y^\n(t>'{U) (9) 



where L^^^ is the luminosity limit of the survey at red- 
shift z. The function (/''(L) was then varied to map out 
the distribution of likelihoods as a function of logL*, a, 
and the mean density of galaxies brighter than given line 
luminosity. Note that this style of analysis is similar to 
the classical method of Sa ndage et al . (1979), with the 
double integral ter m providing the normalization for (j) 
prorvet alj[2003l) . 

Figure ini shows the likelihood contours for the observed 
luminosity function. The left-hand panels marginalizes 
over the integral of the luminosity function (which serves 





Figure 6. Likelihood contours of our lSchechterl I I1976I1 fits to the 
observed [O II] A3727 emission-line luminosity functions for our 
four redshift bins. The contours are drawn at 0.5 cr intervals from 
0.5 cr to 3(T. The left-hand panels marginalize over normalization 
and show probability versus logL* and a\ the right-hand pan- 
els marginalize over a, and plot likelihoods versus logL* and the 
integral of the Schechter function. (The latter parameter is less 
dependent on a and L* .) Note the increase in L* with redshift 
and the poorly constrained measurement of a. 

as a proxy for (jf ) and shows the probabilities as a func- 
tion of a and logL*; the right-hand side marginalizes 
over a. The first feature to notice is that the faint- 
end slope of our function is poorly defined; this is a 
natural consequence of the limited number of objects 
in our sample. Literature estimates of a vary substan- 
tially, from a = -1 2 JGallcgo et al. 2002]) to a = -1.6 
(jSuUivan et al.l[2000h . but our data imply slopes that are 
even shallower. However, this result is quite sensitive to 
the incompleteness corrections which have been applied 
to the measurements near the ^ 5 u detection limits. For 
the discussion that follows, we fix a = —1.2, while noting 
that shallower (or steeper) slopes cannot be excluded. 

Table [1] lists of best-fitting values of L* and (jf under 
the assumption that a — —1.2. These fits are also dis- 
played in the left-hand panel of Figure [SJ In the two 
lowest redshift bins [z < 0.325) where our survey vol- 
ume is relatively small, the bright-end of the luminosity 
function is sparsely populated, and the exact location of 
L* is defined as much by our assigned value of a and 
by the absence of luminous objects as the presence of 
L > L* galaxies. Nevertheless, the exponential cutoff in 
the lowest redshifts bins is reasonably well-defined with 
logL* = 41.07±0.17atz < 0.2 and logL* = 41.29±0.11 
at 0.2 < z < 0.325. These values are significantly fainter 
than those measured in the higher redshift bins, where 
our survey volume is larger, and the knees of the lumi- 
nosity functions are better populated. Overall, the data 
imply an increase in L* of ~ 1.5 mag per unit redshift 
over the redshift limits of our survey. 

Figure [7] shows the results of our maximum- likelihood 
fits after de-reddening the [O II] A3727 line-luminosity 
of each galaxy using our SED-based extinction values. 
Predictably, the evolution of L* with redshift is stronger 
than for the raw luminosity function. As pointed out 
by a number of authors (e.g., ,Wang fc Heckman 19961: 





Figure 7. Likelihood contours of our Schechteri (1976;) fits to the 
de-reddened [O II] A3727 emission-line luminosity functions for our 
four redshift bins. No correction for galaxies extinguished below 
the survey's completeness limit has been applied. The contours 
are drawn at 0.5 a intervals from 0.5 cr to 3 it. The left-hand panels 
marginalize over normalization and show probability versus log L* 
and a; the right-hand panels marginalize over a, and plot like- 
lihoods versus logL* and the integral of the Schechter function. 
Because extinction correlates with luminosity, the evolution of L* 
with redshift is more dramatic than for the observed luminosity 
function. 



iHopkins et all 120011: iLv et al.|[20l2l ). the amount of ex- 
tinction affecting a galaxy's emission lines strongly corre- 
lates with its star formation rate, with more active galax- 
ies possessing more internal extinction. Consequently, as 
the knee of the observed luminosity function shifts to- 
wards brighter values, the amount of internal reddening 
increases, leading to a larger shift in the de- reddened 
value of L*. At 2 < 0.2, the difference between the 
best-fitting observed and de-reddened values of L* is 
AlogL* = 0.4 ± 0.2 dex; by z ~ 0.5, this offset be- 
comes at least AlogL* = 0.7 ±0.2 dex. The brightening 
of L* with redshift increases the importance of internal 
reddening, and explains why our de-reddened values for 
L* are brighter than those of iLv et aD (|2007| ). who used 
a mean reddening correction for all the galaxies in their 
z ^ \ survey. 

We do note that our de-reddened luminosity function 
displayed in the right-hand panel of Figure [5] is almost 
certainly an underestimate. While one can count the 
observed number of galaxies in some interval about a de- 
reddened luminosity Lq) there is always the possibility 
that some additional galaxies which belong in the inter- 
val are extincted below the flux limit of the survey. Obvi- 
ously, the problem becomes more serious at fainter mag- 
nitudes, but at higher redshifts, even relatively bright 
galaxies may be missing from the data. 

To examine the possible importance of this effect, we 
can use the MPA-JHU spectrum measurements of galax- 
ies in the SDSS Data Release 70. This catalog contains 
the line fluxes and Balmer-line extinction estimates for 
several hundred thousand galaxies in the nearby uni- 
verse, allowing us to examine the distribution of [O II] 
extinctions as a function of intrinsic [O II] luminosity. 

® http://www.mpa-garching.mpg.de/SDSS/DR7/ 



These data illustrate a smooth trend in both the mean 
extinction of a population and the dispersion about this 
mean. For example, for systems with [O II] luminosities 
of ~ 10^^ ergs s~^, the mean extinction is E{B — V) ~ 0.2 
with a dispersion of <Je(b-v) ^ 0-13 about this mean. 
By L{[0 II]) ~ 10''3 ergs s'^ however, {E{B - V)) has 
increased to ~ 0.75 and the dispersion has increased to 
<^E{B-v) ~ 0.36. If this relation holds through z ~ 0.5, 
then, at any luminosity and distance, we could, at least 
in theory, compute the fraction of galaxies that should be 
extinguished below our completeness limit, and include 
that as an additional factor in the computation of the 
modified luminosity function 4)'{L). 

Unfortunately, the current sample of galaxies is not ex- 
tensive enough to apply this analysis technique or test 
the validity of our assumptions about extinction. To 
study galaxies beyond the knee of the z > 0.4 [O II] 
A3727 luminosity function, one must reach intrinsic [O II] 
luminosities of ~ 10*^ ergs s~^. Such objects have a mean 
internal extinction of E{B — V)^ 0.4, but with a disper- 
sion of <te{b-v) ^ 0.15 about this mean. Consequently, 
in order to keep the incompleteness corrections below 
~ 50%, one must reach observed [O II] luminosities of 
~ 10^° ergs s~^, and exclude detections within ^ 1 dex 
of the frame limit. At redshifts z > 0.3, our current data 
do not have a sufficient number of galaxies to satisfy this 
criteria. The full HETDEX survey will certainly be able 
to investigate this question. 



7. 



EVOLUTION OF THE STAR FORMATION RATE 
DENSITY 



For individual galaxies, the relationship between [O II] 
luminosity and star formation rate is complicated, since 
it depends on variables such as extinction, metallic- 
ity, and excitation. However, relatively robust calibra- 
tions for galaxy ensembles do exist in the literature 
(iJansen et al.ll2001l:lKewlev et al.ll2004l: iMoustakas et al.l 
l2006f ). and these can be used to translate our [O II] 
luminosity functions into star formation rate densities. 
As Figure [4] illustrates, the lKewlev et al.l (|2004[ ) relation 
(without its second-order metallicity correction) gives 
star formation rates that are reasonably consistent with 
measurements in the UV. This is the conversion we adopt 
for our study. 

Table [l] hsts two estimates for the [O II] A3727 star 
formation rate density: one derived by extrapolating a 
faint-end slope of a = —1.2 to infinity and applying a 
sing le mean extinction of ^.-^727 = 1.88 rn ag to all galax- 
ies (|Hopkinsll200l iTakahashi et al.ll2007[ ). and a second 
using the same faint-end extrapolation, but with our 
individual extinction estimates. (Our luminosity func- 
tions are sufficiently deep that this extrapolation only 
increases the observed SFRD by ^ 10%.) In general, the 
SFRDs derived from the individual reddening values are 
'^0.3 dex smaller than those inferred using a mean ex- 
tinction. This should not be too surprising: while the 
latter method more accurately corrects for extinction at 
the bright end of the luminosity function, it most likely 
underestimates the function's faint-end slope, due to the 
contribution of galaxies whose internal extinction causes 
them to fall out of our sample. Conversely, the applica- 
tion of a mean extinction likely overestimates the faint- 
end slope, since the E{B — V) correlates directly with 
line luminosity. A steepening of a by just ~ —0.4 would 
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Figure 8. Extinction-corrected measurements of the cosmic star 
formation rate density based on the emission- lines of Ha (in red), 
[O III] A5007 (in green) and [O II] A3727 (in blue). The circles show 
the data from this paper. Included in the p lot are measurem ents 
fro m th e literature using the calibrations of Kennicutt (1998) for 
Ha.lKewley et al. (2001) for [O II] A3727, and .Lv ct al.. 12007.) for 
[O III] A5007. 

move the two values into agreement. 

Figure E] compares our individually de-reddened 
SFRDs to other emission-line based SFRDs, usi ng the 
calibrations of Kennicutt (1998) for Ha, Kcwlc vet al.l 
(jl004) for [O II] A3727, and lLvet al.l (l2007f l for [O III] 
A5007. Clearly, the diagram has quite a bit of scat- 
ter, in part due to the different assumptions investiga- 
tors make about internal extinction. Nevertheless, the 
well-known correlation between star formation rate den- 
sity and redshift is easily seen, and our data agree with 
those of most other measurements of instantaneous star 
formation rate. At z '^ 0.5, our value of psfrd = 
3At\i X 10 



Mq yr JVIpc •^ is ^ 4 times greater than 
that seen in the local (z < 0.2) universe, yet a factor 
of ~ 2 below that inferred from [O II] observations at 
z ^ 1. Our data populate the gap between the [O II] 
measurements in the local and high-redshift universe. 

8. CONCLUSIONS 

We have used the sample of 284 [O II] -selected galax- 
ies from the HETDEX pilot survey to explore evolution 
in star formation over the last ~ 5 Gyr of cosmic time. 
By analyzing the distribution of HETDEX line fluxes, 
we find that, over the redshift range < z < 0.56, 
the observed [O II] luminosity function exhibits signif- 
icant evolution, with log L* fading by ^ 0.4 dex between 
z ~ 0.5 and z ~ 0. However, this shift tells only part 
of the story. Because the HETDEX pilot survey was 
conducted in areas of the sky with large amounts of an- 
cillary data, we used each galaxy's spectral energy dis- 
tribution to estimate its internal reddening. By adopt- 
ing these stellar E{B ~V) measurem ents, an d assuming 
E{B - y)stars = 0A4:E{B - y)gas (jCalzettii .2001) . we 
were able to correct each galaxy's observed [O II] flux 
for internal extinction, and thereby obtain an estimate 
of its intrinsic [O II] luminosity. The resultant extinction- 
corrected luminosity function exhibits stronger evolution 
in log L* than the observed function, with this charac- 
teristic luminosity fading by ~ 1.3 mag between z '^ 0.5 
and today. The difference between the observed and 
de-reddened functions is consistent with the expecta- 
tion that galaxies with higher star formation rates have 
greater internal extinction. 



A similar evolutionary trend is seen in the strength of 
the [O II] line relative to the continuum. While the distri- 
bution of rest-frame [O II] equivalent width is insensitive 
to line-luminosity, it does depend strongly on redshift. 
Specifically, there are many more high-equivalent width 
[O II] emitters at z ~ 0.5 than there are today. If we 
exclude the weak line emitters (i.e., those objects with 
rest-frame equivalent widths less than 5 A) and fit an 
exponential to the remainder of the distribution, then we 
find that the scale length of rest-frame equivalent widths 
decreases smoothly, from Wq = 22 ± 3 A at z > 0.45 to 
Wo — 8±2 A at z < 0.2. This reddening- independent 
measure of star formation confirms that the relative in- 
tensity of galactic starbursts has been decreasing over 
the last ^ 5 Gyr. 

By integrating our extinction-corrected [O II] luminos- 
ity functions and u sing the [O II] SFR calibration of 
IKewlev et al.l (|2004l ). we derive the star formation rate 
density in four redshift bins. We find that the SFRD 
decreases linearly with time, changing by a factor of ~ 4 
between z '^ 0.5 and today. Our data cover the gap be- 
tween the [O II] observations of the local universe and 
those at z ~ 1, and are in excellent agreement with mea- 
surements based on the Ha and [O HI] A5007 line. There 
remains a substantial amount of scatter in the SFRD di- 
agram, but this is due in large part to differences in tech- 
niques and reddening assumptions. Our analysis, which 
employs individual SED-based reddening estimates and 
spans the redshift range < z < 0.5, removes some of 
these problems. 

The sample presented here highlights the pow er of the 
upcoming HETDEX survey (|Hill et a l." '2008a') to open 
up the emission-line universe. HETDEX will map out 
over 300 deg^ of sky with a blind, integral-field spectro- 
scopic survey. While the main goal of the project is to 
measure the power spectrum of ~ 800, 000 Lya emitting 
galaxies between 1.9 < z < 3.5 and measure the evolu- 
tion of the dark energy equation of state, the survey will 
also identify ^ 10® z < 0.5 [O II] emitters. This unprece- 
dented large sample of emission-line galaxies will cover 
the entire range of galactic environments and allow us to 
trace the evolution of star formation in field galaxies and 
clusters over the last ^ 5 Gyr of cosmic time. 
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Table 1 

Best-Fit Schcchtcr Function Parameters^ 



Observed Extinction-Corrected 

Wo log-L* log </>(> 40.5) logpsFRD logL* log </>(> 40.5) logpsFRD 

Nahj (A) (ergs s^l) (Mpc-3) (^Mq yr~l Mpc-^) (ergs s^l) (Mpc^^) (Mq yr^^ Mpc^^) 

'-1:11 -2.30in? -1-63 ±0.11 41.48l»:?J -2.02l°;E;i 

0.200-0.325 70 11.5 ± 1.6 41.291!];" -'^■'^'^-om -1.40 ±0.06 il.dG'ho^^ -'^■^'^-o.ot -1.67 ±0.08 



0.000-0.200 39 8.0 ±1.6 41.071016 -2.30Ioii -1.63 ±0.11 41.48Ioi3 -2.02I'^;;g -1.99 ±0.11 

t»ll -2.12lHt -1.40 ±0.06 41.960-12,, -1.921^5 

0.325-0.450 89 16.6 ± 2.4 41.50t„;?o ~2.07to:E55 -1.29 ±0.05 42.160-J°ii -'^■^^to.os -1.60 ±0.06 

0.450-0.560 76 21.5 ± 3.3 41.68t„;"2 -2.07to;°| -1.24 ±0.06 42.36'ljii2 -2-02lo:o7 -1.56 ±0.07 

^ Parameters assume the local value of a — —1.2 



